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Abstract 

We consider space-times with two isometries which represent gravitational waves 
with distinct wavefronts which propagate into exact Friedmann-Robertson- Walker 
(FRW) universes. The geometry of possible wavefronts is analysed in detail in all 
three types of FRW models. In the spatially flat and open universes, the wavefronts 
can be planar or cylindrical; in the closed case they are toroidal. Exact solutions are 
given which describe gravitational waves propagating into the FRW universes with a 
fluid with a stiff equation of state. It is shown that the plane-fronted waves may include 
impulsive or step (shock) components, while the cylindrical waves in the spatially flat 
and open universes and the toroidal waves in closed universes may contain steps. In 
general, wavefronts may exist which have an arbitrary finite degree of smoothness. 
In all cases, the waves are backscattered. The head-on collision of such waves is also 
briefly mentioned. 
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1 Introduction 



Most work on gravitational waves in cosmological backgrounds has been based on first 
order approximation methods. However, if gravitational radiation played a significant role 
in the universe, then the waves would cause the symmetries of the standard homogeneous 
and isotropic models to be violated. In order to treat general, large perturbations one 
has to turn to a numerical approach. However, when space-time is assumed to have two 
spacelike Killing vectors, exact solutions containing gravitational waves or solitons can be 
constructed, at least for vacuum or for scalar or electromagnetic fields or for a stiff perfect 
fluid in which the pressure is equal to the density. Some of these solutions have been 
reviewed for example by Carmeli et al 0], by Adams et al |^], and recently by Verdaguer 

In previous work on gravitational waves in cosmological backgrounds, the exact solu- 
tions considered have been interpreted as perturbations extending over whole universes. 
As such, it has not been possible to give an analysis of the geometry of the wave surfaces. 
In this paper we present the new class of solutions representing gravitational waves prop- 
agating into homogeneous, isotropic universes. The space-time ahead of a gravitational 
wave is taken to be an exact Friedmann-Robertson- Walker (FRW) universe. Since this 
has spatial sections of constant curvature, the geometrical properties of the wavefront can 
be explicitly determined. For these solutions, we can also analyse the character of the 
wavefront: whether it contains a ^-function impulse, or a shock (step) wavefront, or is 
smooth in the sense that the first n derivatives of the Weyl tensor are continuous across 
the front. 

We consider space-times which admit two commuting, hypersurface orthogonal space- 
like Killing vectors. These can be described by the line element 

ds 2 = e - M (d V 2 - dfi 2 ) - e- u (e v dx 2 + e" v ' dy\ (1.1) 

where the metric functions U, V and M depend on the coordinates rj and fi only. In 
fact, all the FRW models can also be written in this form for a perfect fluid with any 
equation of state, with or without a cosmological constant. Thus, certain exact solutions 
described by this line element can be considered to be (large) perturbations of these FRW 
models, provided they reduce to particular models in appropriate limits. (Other classes 
of inhomogeneous universes which can reproduce the FRW models in certain limits have 
recently been reviewed by Krasihski 

It may immediately be noticed that, in the form of the line element (|1.1|), two families 
of null hypersurfaces, ry ±/i = const., are singled out in a natural way. These surfaces may 
be regarded as wave surfaces for both approximate and exact gravitational waves. The 
geometry of these wavefronts in the FRW universes of all three types will be analysed in 
detail in this paper. These null hypersurfaces may also be considered as suitable wavefronts 
for some other radiative fields. 

In a series of recent papers some families of exact solutions have been outlined 

which describe the propagation (and collision) of gravitational waves into FRW back- 
grounds in the above form. In these solutions, the space-time is exactly FRW with stiff 



2 



fluid ahead of the waves. Such waves therefore have well-defined wavefronts. It was pre- 
viously noted that these waves are necessarily backscattered, and that their interactions 
following head-on collisions do not induce additional singularities into the space-time. 
However, complete metrics were not given and the structure of the wavefront was not 
analysed. A method for obtaining complete solutions was later described in ||. 

It is the purpose of this paper to complete the derivation of the solutions representing 
waves with distinct wavefronts propagating into FRW backgrounds, and to summarise 
them from a unified point of view for each possible type of wavefront and all three types 
of FRW models. We will describe the waves by considering the most general form of the 
solution of the main field equation for a distinct wavefront. Using the method outlined 
in H, we will complete the integration of the subsidiary equations, thus determining the 
whole metric. We will also analyse the structure of the wavefront in each case. It will be 
shown that these may include impulsive components in some cases, but may contain step 
(or shock) waves in all cases. Waves in which the Weyl tensor has n continuous derivatives 
are also discussed. 

We will first describe the geometrical character of possible wavefronts in each FRW 
model in order to treat the properties of all the above solutions from a unified point of view. 
In Sections 2-4, we will consider the spatially flat, open and closed FRW models in turn. 
It will be shown that, in the spatially flat and open FRW universes, the gravitational 
waves can have plane wavefronts. In these universes, the gravitational wavefronts can 
alternatively be cylindrical. In the case of the closed FRW universe, the wavefronts are 
toroidal. In the Appendix, we describe in detail how a 2-torus, forming a wavefront, can 
be constructed in the closed model. The character of the wavefronts reflects only the 
symmetries assumed. This approach is thus relevant also in any study of more general 
waves and cosmological models than those treated here. It may be of interest to analyse 
even small (linear) perturbations of FRW universes with these symmetries, with other 
equations of state and, possibly, with a non-vanishing cosmological constant. 

After summarizing the properties of the FRW universes with a stiff fluid in Section 5, 
we study in Section 6 the general family of solutions which describe waves with distinct 
wavefronts propagating into these backgrounds. Then, in Sections 7-11, we analyse re- 
spectively the separate cases of plane and cylindrical waves in spatially flat and open 
universes, and toroidal waves in a closed universe. 

As well as analysing the character of the wavefront, we also find that, for the waves 
with plane wavefronts the only space-time singularity is that describing the big bang. 
However, our cylindrical and toroidal waves will, in general, propagate away from (and 
be backscattered towards) a singular axis, which can be considered as both "source" and 
"absorber" of the waves. Other remarks appropriate by way of conclusion are collected in 
Section 12. 
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2 Wavefronts in spatially flat FRW universes 



In the spatially flat case (with curvature index k = 0), the Robertson- Walker line element 
may be written in either of the forms 

ds 2 = dt 2 - R 2 {t)(dx 2 + dy 2 + dz 2 ), (2.1) 
ds 2 = dt 2 - R 2 {t){d P 2 + p 2 d(f) 2 + dz 2 ), (2.2) 

where R(t) is determined by Friedmann's equation. In order to consider null hyper surf aces, 
it is convenient to introduce the "conformal time" coordinate r] such that 

dt = R(r])d7]. (2.3) 



The line elements (2.1) and ( |2.2| ) can then be written in the form ( |1.1|) : 



ds 2 = R 2 {rf){dr] 2 - dz 2 ) - R 2 {r,){dx 2 + dy 2 ), (2.4) 
ds 2 = R 2 ( V )(dr ] 2 -dp 2 )-R 2 (r ] )(p 2 d^ 2 + dz 2 ). (2.5) 

In (|2.4|), surfaces on which rj = const., z = const, are planes, and the null hypersurfaces 
rj — z = const, represent plane wavefronts. Clearly, plane wave surfaces oriented in any 
direction can immediately be constructed. Alternatively, in the form ( |2.5| ), surfaces on 
which i] = const., p = const, are cylinders, and the null hypersurfaces rj — p = const, 
represent expanding cylindrical wavefronts. 

3 Wavefronts in FRW open universes 

The standard form of the metric of an open FRW model with negative spatial curvature 
(k = —1) in coordinates x £ [0,+oo), 8 £ [0,7r), </> G [0, 2tt) reads 

ds 2 = dt 2 - R 2 {t)[d X 2 + sinh 2 X (d6 2 + sin 2 6 dcf 2 )}. (3.1) 

As is well known | |10[ |) the spatial 3-geometry on a hypersurface t = to = const., 

dt 2 = R 2 [d X 2 + sinh 2 X {d9 2 + sin 2 6 d<j> 2 )) , (3.2) 

where Rq = R(to), can be represented as 3-dimensional hyperboloid 

V 2 - X 2 -Y 2 - Z 2 = R 2 , (3.3) 

embedded in a 4-dimensional Minkowski space with Lorentzian coordinates V, X, Y, Z, 
where 

V = .Rocoshx, 

Z = i?o s i n hx cos $; 

X = i^osinhxsin^cos^, 

Y = -Roshih^sin^sinc/;. (3-4) 
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In order to find 2-planes in the open FRW model, we introduce alternative coordinates 
in which the 3-metric is conformally flat. By putting 

x = Ro—, y = Ro—, 

z z 
Z = ^(l-x 2 -y 2 -z 2 ), (3.5) 

V = ^(l + x 2 + y 2 + z 2 ), 



where x, y € (—00, 00), z S (0, +00), we find Q3,3j ) to be satisfied. These coordinates sweep 
out the whole "positive" (V > 0) sheet of the hyperboloid (|3.3j). The 3-metric now takes 
the form 

dl 2 = \(dx 2 + dy 2 + dz 2 ). (3.6) 

z 

Clearly, putting z = const, we obtain 2-surfaces with intrinsically flat geometry. Both the 
3-metric ( |3.6| ) and the whole space-time metric admit, on these 2-surfaces, not only the 
translational Killing vectors d x and d y , but also the rotational Killing vector yd x — xd y . 
The proper lengths along the whole orbits of both d x and d y are infinite. Thus, these 
2-surfaces on which z = const, can be considered to be 2-planes. 

In figure 1 we illustrate these 2-planes in the embedding diagram of the section Y = 0, 
or (f> = and tt, through the 3-hyperboloid ( |3.3| ) in a 3-dimensional Minkowski space with 
coordinates V, X, Z. Since z = zq = const, implies V + Z = Rq/zq = const., the 2-planes 
are formed by cutting the 3-hyperboloid ( p.3[ ) along "null hyperplanes" V + Z = const, 
in the Minkowski space. In this way we obtain the 2-planes described in coordinates 
V,X,Y,Z by 

-Z = i ZQ (X 2 + Y 2 ) + Z , 
Zq = \zqRq 2 — zq^ 1 = const., (3-7) 
V = Rqzq^-Z. 

For a fixed zo, this is a 2-paraboloid which in figure 1 reduces to a parabola, the Y 
dimension having been suppressed. (Figure 1 remains the same after replacing X by Y .) 
On the other hand, by considering sections of Minkowski space on which Z = const., 
equations (3J3) and ( |3~7^ ) for fixed z = zq imply X 2 + Y 2 = const., or x 2 + y 2 = const. 
With these cuts of the 3-hyperboloid, the infinity of the 2-planes cannot be seen. Indeed, 
in terms of the original coordinates %■> t ne condition V + Z = RqZq -1 reads cosh \ + 
sinhxcos# = zq~ 1 . In order to consider \ ~^ °°i we have to let 6 — > tt, as expected for 
paraboloids with axes in that direction. In choosing coordinates according to (|3.5| ), we 
"singled out" the coordinate z; since the 3-hyperboloid is isotropic and homogeneous, we 
can construct other sets of 2-planes by simple transformations. 

Cylindrical 2-surfaces in the open FRW models can easily be found by introducing 
coordinates p G [0, +00), y S (— 00, +00), (p G [0, 2ir) such that 

V = Ro cosh p cosh y, 
Z = Rq coshpsinhy, 
X = Rq sinh p cos </>, 

Y = Rq sinh p sin (f>. (3.8) 
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The 3-metric on the 3-hyperboloid then takes the form 



d£ 2 = R 2 [dp 2 + sinh 2 pdcp 2 + cosh 2 pdy 2 }. 



(3.9) 



Putting p = po = const., we obtain cylindrical 2-surfaces. The proper length along the 
orbits of the Killing vector cL is finite, being 2ttRq sinh p, but it is infinite along the orbits 
of d y . Although the local geometry on the cylinder is intrinsically flat, there exists no 
other Killing vector in this 2-space as occurs in the case of the 2-plane. 

In figure 2, a cylindrical 2-surface p = po is illustrated in the same embedding diagram 
as in figure 1. From (3.8), we now find 

X 2 + Y 2 = (i? sinhp ) 2 = const., 



V 2 - Z 2 



(i?o cosn Po) 2 = const. 



(3.10) 



Since the Y-dimension is suppressed ((f) = and %) in figure 2, we obtain X = ±i?o smh po 
and V 2 — Z 2 = (Rocoshpo) 2 . These are two hyperbolae, representing two generators of 
the cylindrical 2-surface p = po. The axis of the cylinder is given by p = 0, i.e. 9 = or 
7r, or X = Y = 0. Considering the section of the 3-hyperboloid ( |Q| ) on which Z = 0, the 
2-cylinder would be represented by the circle X 2 + Y 2 = (Rq sinh po) 2 , V = Rq cosh po, as 
follows from ( |3.10| ). Clearly, for po small, the 2-cylinder would approximate to an ordinary 
cylinder in E 5 . 

In terms of the two new sets of coordinates ( |3.5D and ( |3.8| ), the line element fl3.1| ) takes 
the following forms 



ds 2 
ds 2 



dt 2 - R 2 {t) 



-^(dx 2 + dy 2 + dz 2 



dt 



R 2 (t) dp 2 + sinh 2 p d(f> 2 + cosh 2 p dy' 



Introducing the new coordinate 



the line element (3.11) becomes 



p = logz, 



ds z = di l - R\t) dp 1 + e~ 2il {dx l + dy 1 ) 



(3.11) 
(3.12) 

(3.13) 
(3.14) 



Finally, going over to the "conformal time" coordinate r/ in both cases, we arrive at the 
metric 

ds 2 = R 2 (r])(dri 2 - dp 2 ) - R 2 (r])e- 2fl (dx 2 + dy 2 ). 



Alternatively, the line element ( p,12| ) can be written as 

ds 2 = R 2 (r]){dri 2 - dp 2 ) 
Both metrics ( |3.15| ) and 3.16) are now expressed in the form (1.1). 



R 2 (rj) sinh 2p( t&nh pd(p 2 + coth pdy 2 ) 



(3.15) 
(3.16) 
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4 Wavefronts in FRW closed universes 



We start with the standard form of the metric for a closed FRW model having positive 
spatial curvature (k = 1) 

ds 2 = dt 2 - R 2 (t)[d X 2 + sin 2 X {d0 2 + sin 2 6 # 2 )], (4.1) 

where % £ [O? 71 "]) £ [0,tt], ^ £ [0, 27r]. It is well known that the spatial 3- geometry on a 
hypersurface of homogeneity t = to = const., 

d£ 2 = R^dx 2 + sin 2 X {d6 2 + sin 2 9 d<p 2 )}, (4.2) 

where Rq = R(to), can be visualised as a 3-dimensional sphere 

W 2 + X 2 + Y 2 + Z 2 = R 2 ,, (4.3) 

embedded in a 4-dimensional Euclidean space with cartesian coordinates W, X, Y, Z. By 
setting 

W = Rq cos 
Z = i?osinxcos#, 
X = Rq sin x sin 9 cos </>, 

V = Rq sin x sin sin </>, (4.4) 

Eq. ( p~3"D is satisfied. 

In the Appendix it is demonstrated in detail that, with the symmetries assumed, a 
typical wave surface in a closed FRW model is a 2-torus in the 3-sphere. This is best 
described by parametrizing the whole 3-sphere by coordinates Q € [0, vr/2], a E [0, 2tt], 
5 € [0, 2ir] such that (cf. Eq. flAjf ) in the Appendix) 

W = Rq cos C cos a, 

Z = Rq cos C sin cr, 

X = Rq sin (" cos <5, 

Y = i? sinCsin(5. (4.5) 

The 2-torus is specified by the choice of the parameter k = cos £. (For an illustration, see 
Figure 8 in the Appendix.) 



In terms of the new coordinates, the line element (4.1) takes the form 

ds 2 = dt 2 - R 2 (t)[d( 2 + sin 2 (d5 2 + cos 2 (da 2 }, (4.6) 
or, using the confer mal time r], 

ds 2 = R 2 (7])(drj 2 - d( 2 ) - ii? 2 (r?)sin2C(tanC^ 2 + cot (da 2 ), (4.7) 



which is clearly in the form (1.1) 
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5 FRW universes with a stiff fluid 



We now consider the case of a "stiff" perfect fluid in which the pressure is equal to the 
density. This equation of state may be considered as a "modern version of aether" (see |Q]) 
which enables mathematical relativists to generate a variety of new solutions. However, it 
may also have played an important role in the period of high densities after the big bang 



1 11], [12], or even in the interior of cold condensed objects. Its role as a limiting equation 
of state of neutron matter was studied thoroughly in the full field theory by Walecka 
Jk|] . Difficult and important problems arise, of course, when studying gravitational waves 
which propagate in dust or in fluids with a more realistic equation of state. However, in 
the case of a stiff fluid, families of exact radiative solutions can be obtained which can be 
considered to be perturbations of the FRW models with a stiff fluid. In this section, we 
review the unperturbed background models. 

The equation of state of the fluid is taken to be p = p. (We use tildes to denote 
quantities referring to the fluid to avoid confusion with the same symbols used elsewhere.) 
With this, the equation of energy conservation {pR?),t = —p{R?),t implies that 

pR 6 = p R§ = const. = 70 = ^ 7 2 , (5.1) 

where the dimension of the constant 7 is (length) 2 . The standard Friedmann equation, 
R, t 2 = -k + 7 2 iT 4 , k = 0, ±1, can then easily be integrated after introducing the confor- 
mal time coordinate 77 in accordance with ( |2.3| ), The solutions read ][14|], |7| 

R = [7sinh2?7]5, if h = — 1, 

R = [ 7 2?7] I, if k = 0, (5.2) 

R = [7 sin 2rj\ if k = +1. 



In Q the physical time t is expressed explicitly in terms of 77 by means of elliptic 
integrals in the cases k = ±1. For the spatially flat case one gets simply 



i?=(3 7 t)3, (5.3) 

so that the flat FRW universe with a stiff fluid expands slower than the corresponding 
universe with dust for which R ~ £3 (see e.g. fl(|). For the open universe, one finds 
R = t at large times, i.e. exactly the same asymptote as the standard open Friedmann 
universes with dust or radiation. All these open universes thus approach the flat Milne 
universe, for which R = t exactly. 

The "gravitational role" of the pressure in slowing down the expansion in the spatially 
flat case can also be seen in the "speeding up" of the recollapse in the case of the closed 
FRW models. The conformal time 77 is also often called the "arc parameter" since it 
determines the radius of arc distance on the 3-sphere covered by a photon travelling since 
the start of the expansion [1C]. In the case of the closed universe with dust, the range of 77 



is 2tt — a photon makes "one trip" around the whole universe before its final collapse. In 
the model with radiation, in which p = p/3, the range of 77 is tt — a photon only gets as far 
as the antipodal point of the universe. In the universe with a stiff fluid, the last equation 
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in ( |5,2| ) shows that the range of rj is tt/2 — a photon succeeds in traveling only one quarter 
of the distance around the whole universe (say, from the pole to the equator). This will 
be relevant in Section 11 where we shall consider exact toroidal waves propagating into a 
closed FRW background with a stiff fluid. 



It is easy to write the line element of all three types of FRW models in the form (1.1). 
Since 7 is an arbitrary constant with dimension (length) 2 , we may consider (1.1) in the 



dimensionless form ds 2 /^ 2 . Without loss of generality we may thus put 7 = 1. Comparing 
now (pip, (^), ( gig ), pl6| ) and |L7|) with R(rj) substituted from (jj^) with the general 



form ( |1 . lp , we find the following: 

Case la. In the spatially flat case with the plane 2-surfaces z = const., we have p = z G 
(— 00, +00), (x,y) G IR 2 , 77 G [0,oo), and 

e~ u = 2rj, e' M = 2rj, V = 0. (5.4) 

Case lb. In the spatially flat case with the cylindrical 2-surfaces p = const., p = p G 
[0, 00) , x = <f> G [0, 2-7t) , y = z £ (— 00, +00) , 77 G [0, 00) , and 

e" 17 = 2r/p, e~ M = 2tj, e v = p. (5.5) 

Case 2a. In the open model with the plane 2-surfaces z = const, or p = const., we have 
p G (—00, +00), (x,y) G IR 2 , T] G [0,oo), and 

e~ u = sinh 2r? e" 2/ \ e~ M = sinh 2rj, V = 0. (5.6) 

Case 2b. In the open model with the cylindrical 2-surfaces p = const., p = p G [0, 00), 
x = (J) G [0, 27r), y = z G (—00, +00), 77 G [0, 00), and 

e~ u = sinh 2t? sinh 2p, e~ M = sinh27/, e^ = tanhp. (5-7) 

Case 3. In the closed model with the toroidal 2-surfaces £ = const., we have p = £ G 
[0,vr/2], x = <5 G [0,2tt), y = a G [0,2tt), 7/ G [0,tt/2], and 

e~ u = sin 277 sin 2(, e~ M = sin 2r/, e y = tan (. (5.8) 



It is well known 15] that, in a space-time with the metric (|1.1]), a stiff perfect fluid 
can be associated with a scalar potential 17(77, p) such that the density and 4-velocity of 
the fluid are given by 



(5.9) 



16tt/3 = e M (a 2 - <r 2 ), u a = a^/{a^a^) 112 , 
and, as a consequence of Einstein's equations, the fluid potential a satisfies 

a m - U v a v - + Un&p = 0. (5.10) 
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For the above FRW universes with a stiff fluid, the expressions for a are found to be 

o- = 3 1/2 log7? for fc = 0, (5.11) 

6r = 3 1 / 2 logtanh»7 for k = -1, (5.12) 

a = 3 1/2 logtanr/ for k = +1. (5.13) 

In the spatially flat and open cases, it should be noted that these expressions for a apply 
both in the plane and cylindrical cases. Clearly, all coordinate systems above are comoving 
with the fluid. 

Using these background metric functions, we now turn to the construction of the waves 
which propagate into these backgrounds. 



6 Waves in FRW universes with a stiff fluid 



We return to the general line-element ( |1.1[ ). In the case of a stiff perfect fluid, and vanishing 
cosmological constant, Einstein's field equations imply that e~ u satisfies the wave equation 

(e- u ) m - (e- u )^ = 0, (6.1) 

and V satisfies the linear equation 

V m - U V V V - + = 0. (6.2) 

It may immediately be noted that, in general, Eq. ( |6.1| ) can be integrated to give 

e~ U = f(v-ti +9(V + »), (6-3) 

where / and g are arbitrary functions. In many cases, it is convenient to adopt / and g 
as coordinates, at least in some region of space-time in which they are not constant. The 
arbitrariness thus retained can be used to describe families of exact solutions. It can be 
seen, however, that singularities of some type will occur when / + g = 0. In addition, in a 
closed universe, / and g cannot be monotonic functions over the entire space-time, and it 
is then more convenient to revert to coordinates that are more directly related to r] and \x. 

Taking / and g as coordinates, equations ( |6.2| ) and ( 5.10 ) take the form of Euler- 
Poisson-Darboux equations with non-integer coefficients 

(f + 9)V f9 + lV f + lV g = 0, (6.4) 
U + g)d fg + \d f + \d g = V. (6.5) 

For any solutions of these equations, the remaining metric function M can be found by 
quadratures. Putting 

( / + 5 )l/2 e - ^ 

the remaining field equations become 

S f = -!(/ + 9)(Vf 2 + Sj), S g = -§(/ + g)(V g 2 + a 2 ), (6.7) 
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which are automatically integrable in view of (|6,4|) and (|6.5D . However, these equations are 
also significant in determining conditions on the derivatives of V and a when considering 
junction conditions across possible wavefronts. In addition, care has to be taken on any 
hypersurface on which /' or g' are zero. Using a null tetrad system that is naturally 
adapted to the coordinates (with the null vectors given by £ a = 2- 1 /2 e -M/2 (?? 
and n a = 2 _1//2 e _M//2 (r/, Q +fJ, a ), the non-zero components of the Weyl tensor are given in 
terms of the coordinates / and g by 



1 j2 " M 1 *v g9 + - (/ + g )v g {v 2 + ~a% 

*2 = -\fg'e M ( 1 j^ 2 ~ W,V B - a f a g ^ , (6.8 



*4 = ~\fe M [iVff + V, - (f + g)Vf(Vj 2 + a}] 



In many cases of interest, it is possible to choose / such that / = represents a null 
hypersurface. On such a hypersurface, r] — /i = const. It is then convenient to adopt a 
gauge such that e~ u = f + g takes the same form on either side of / = 0, thus facilitating 
the possible choice of / and g as coordinates. It may then be noticed that, since (6.4) 



is linear in V, we may superpose a background solution of (6.4) on either side of the 



hypersurface with another solution which is non-zero on one side only. In this way, we can 
construct an exact solution which contains a distinct wavefront which is given by the null 
hypersurface / = 0. Of course, it is necessary to satisfy the appropriate (O'Brien-Synge) 
junction conditions |l^] across the wavefront. In this case, these conditions require that 
both V and M be continuous across / = 0. In the cases we are considering in this paper, 
it is found using Eqs. ( |6.7l ) that the continuity of M implies that Vf must also be bounded 
on / = 0. 

In the class of solutions we wish to consider, the background is taken in turns to be 
appropriate forms of the various FRW models with a stiff fluid given by (^.4|)-( |5^8| ). The 
additional component of V can then be interpreted as introducing a purely gravitational 
wave. In addition, the expression for a is taken to be unchanged across the wavefront. This 
effectively eliminates the introduction of acoustic waves in the fluid of the type considered 
by Tabensky and Taub [19|. However, since Eq. ( |6.5D for a is identical to Eq. ( |6.4D for V, 



acoustic waves could also be considered by including the same kind of additional terms in 
a. Of course, this equivalence occurs only for a stiff fluid in which the characteristics for 
both gravitational and acoustic waves are the same. 

As introduced in [0], we now consider the appropriate class of self-similar solutions of 
( |6.4| ) of the form 

V k (f,g) = (f + g) k H k (f 1 l), (6.9) 

where k is an arbitrary real parameter whose range will be determined below, and the 
function H&(C) is a solution of the equation 

(1 - C 2 )H'l + (2k - l)(H' k - k 2 H k = 0. (6.10) 

This equation admits a class of solutions which can be expressed in terms of hypergeometric 



functions F(a,b;c;z) such that the solutions ( |6.9D of (6.4) can be expressed in either of 
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the forms 



(/ + 9)*ff*0^) = c t l^-F^,l + k;l + k;-i), (6.11) 



where c k = fTfc+^y- Both forms of this expression ( 6.11 ) and (|6.12 ) will be useful below, 



r(fc+3/2) • 

depending on the appropriate ranges of / and g. For integer values of k, these func- 
tions take particularly simple forms in which successive functions H k (() can be obtained 
recursively using the relation 

H k (C)= £ flfc-i(C') dC', (6-13) 

with the initial solution flo(C) = cos -1 ( if / > 0, or i?o(C) = cosh -1 £ if / < 0. It may be 
noted that the scaling constants c k have been included in ( |6.11 ) and ( |6.12| ) to ensure that 



no additional scale factor appears in Q6.13 ), and we may generally put H' k {C) = -fffc-i(C)- 



As can be seen from (6.11) or ( 6, 12] ), the significant feature of the solution fl6,9| ) is that 



V = and Vf is bounded on the null hypersurface / = provided k > |. It may thus be 
matched across the null hypersurface / = to the solution V = on the other side. It then 
represents some kind of strong gravitational wave with wavefront / = 0. As mentioned 
above, since ( |6.4D is linear, this solution may be added to some other background solution 
to represent the propagation of a gravitational wave in that background. For FRW stiff 
fluid backgrounds, the condition k > i ensures that both V and M are continuous across 



the null hypersurface / = 0, thus satisfying the O'Brien-Synge junction conditions [18|. 

In a series of recent papers f|-|ll> fig] , |L7]], families of exact solutions have been 
considered which describe the propagation of gravitational waves in various backgrounds, 
together with the head-on collision of such waves and their subsequent interaction. These 
solutions may here be re-formulated using solutions of (|6.4[) in the form 



roc . . 

V(f,g) = V (f,g) + Q(v ~ m) J i/2 <t>{k)U + 9) k H k (f^) dk, (6.14) 



where Vq(/, g) is a solution of (6.4) representing some particular background space-time, 
@(rj — fj,) is the Heaviside step function in which the wavefront / = is assumed to be given 
by Tj — /i = 0, and (f>(k) is an arbitrary function. The function 4>(k) may be regarded as the 
"spectral amplitude" of some arbitrary wave profile. It is subject only to the condition 
that the above integral exists. 

Let us now turn to the specific cases in which the backgrounds are the various FRW 
models with a stiff perfect fluid source. In these cases, the background regions of course 
are conformally flat. However, the wave regions are algebraically general indicating that 
the propagating gravitational wave is necessarily backscattered. 
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7 Plane waves in a spatially flat FRW model 



In terms of the metric Ql.l| ), the spatially flat FRW universe with a stiff fluid is given by 
and ( |5.11| ). Thus e~ u = 2r], and we can put 

f = r)-z, g = rj + z. (7.1) 

We may thus consider a plane wavefront given by r] = z as described in Section 2, with the 
region r/ < z being an exact FRW background into which a gravitational wave propagates. 

Since V = in the background region in this case, the gravitational wave can be 
introduced by putting 



V( V , z) = 9(7/ - z) Hk)v k H k (J ) dk, (7.2) 



1/2 

in which a factor of 2 k has been absorbed into the arbitrary function 4>(k), and H k (z/rj) is 
given by ( |6.12| ) and its extension. Using the method described in || and |T7j, the complete 
integral for the metric function M can be expressed as 

M(r],z) = -log 2r] + M@(r] - z), (7.3) 

where 

M = -\\ ^dn / 0(fc)0(n - fc) fc(n - k)H k H n ^ k + (1 - 2 T )fl- fc _ 1 fl n _ fc _i rffc, 



(7.4) 

in which the argument of H k (z/rf) has been omitted. 

The non-zero components of the Weyl tensor in this case are given by 

*o = -j-fPWrm + Wvz + V zz - %(y v + Vz )% 

*2 = ^[V r ?-V z \ (7.5) 

*4 = ~^[V m ~ 2V v* + V tt - ~(V V - V z f]. 

These are zero in the background region. Near the wavefront where r\ — z is small, their 
behaviour is governed by the lowest non-zero term in k in ( |7.2[) . Taking this lowest term 
in the form 



V( V ,z) = ?±2 k r, k H k (f j )e(r ] -z) 



= °* {r] F (l' I > I + fe ' K 1 ~ f )) e fa ~ z )> ( 7 - 6 ) 

where a k is a constant, it can be seen that the Weyl tensor components behave near the 
wavefront as 



^~- ^±il fa -^efa- a ), (7.7) 



4, °~-4(lFJ ( ''- 3) ' =+1/2e <"- z )' 
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where 8(rj—z) is the Dirac delta function. This immediately identifies ^4 as the component 
representing the propagating wave with wavefront 77 = z, ^0 is the backscattered wave 
component propagating in the opposite spatial direction, while ^2 represents the effective 
interaction between these two waves. It can also be seen that an impulsive gravitational 
wave is included if the lowest term in the expansion for V has k = ^- It has a step (or 
shock) gravitational wave if the lowest term has k = |. In general, the Weyl tensor is C n 
across the wavefront if k m - m = n + | . 

It can be shown that this solution is regular everywhere except at the initial "big 
bang" singularity 77 = 0, where both the waves and the background may be considered to 
originate. 

These solutions were first considered in || for the case in which k takes half integer 
values only. However, the complete solution was not obtained and the possible existence 
of an impulsive component was not recognised. The complete solution for a sum of terms 
with integer values of k only was given in ||. The head-on collision of such waves was 
discussed in ||] using half integer values of k only, but an explicit expression for the 
function M in the interaction region cannot be obtained by the methods considered so far. 



8 Cylindrical waves in a spatially flat FRW model 

The above method can also be used to construct cylindrical gravitational waves in a 
spatially flat FRW universe with a stiff fluid by expressing it using the coordinates defined 



in case lb of Section 5. In this case, the metric functions are given in the form (5.5) and 

-u 



the fluid by ( 5.11 ). Since e 



/ 



2rjp in this case, we can put 
-Uv-P) 2 , g = Uv + pf- 



.1) 



We may now consider a gravitational wave with a cylindrical wavefront given by 77 = p 
as described in Section 2. Accordingly, the region 77 < p is taken to be the exact FRW 
background into which the wave propagates. Such a wave may be considered to originate 
on a line p = at the initial singularity 77 = 0. In this case V = log p in the background 
region, and a gravitational wave in the general form ( |6.14 ) can be introduced by putting 



V(rj,p) = log p + 6(77 - p) 



1/2 



2 V p 



dk, 



in which a factor of 2 k has again been absorbed into the arbitrary function 4>(k). (This 
function, representing a spectral amplitude, should not be confused with the cylindrical 
coordinate being used in this case.) The functions H^{ p-^- ) are most conveniently given 
here by ( |6.11 ). Putting 



V = logp + VO(r] - p), M = - log 277 + M 6(77 - p), 



(8.3) 



it can be shown that the subsidiary equations ( p.7|) can be fully integrated to yield 



M 



\ \ <p(k)v k P k 

Jl/2 



(v 2 -P 2 ) 1 

2»7p k 



Hk-i 



dk 
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n 



dn 



n-l/2 



1/2 



4>{k)(f){n — k) 



k(n - k)H k H n ^ 



4?7 2 p 2 



dk, (8.4) 



in which the argument of H k is taken to be ^ ^ p p ■ 

In this case, the non-zero components of the Weyl tensor in the region r/ > p are given 

by 



*4 



J_ M 

8r? e 



3rj 
2(»7+p) 



VP 
2(»7+p) 



+ ^p 



'/'/ 



2K 



p^p 



(8.5) 



/'/' 



+ 



v-p 



(V v - V p ) 



Vp? + 2^(V V 



Vo 



The term — g^e (V^ — Vp)<5(r/ — p) which would normally appear in the expression for ^4 



It follows that 



has been omitted because, in this case, V v = V p on the wavefront rj = p. 
this class of solutions does not include impulsive cylindrical waves. 

Near the wavefront where 77 — p is small, the character of the gravitational wave is 
determined by the lowest non-zero term in k in (|8.2| ), Taking this lowest term in the form 

v( v ,z) 



^2 2k V k p k H k (z 

Ck 



■ 2 +p 2 

2rjp 



(v- p) 



l+2k 



iv + p) 



.6) 



it can be seen that, near the wavefront, the Weyl tensor components behave as 



^4 ~ 



-\{2k + 1)(4* + 1) a k rf 2 (rj - p) 2k ~ x @(rj - p), 
-±(2k + l)a k v- 3 (v-p) 2k &(r)-p), 



(8.7) 



pf k+l Q{ri-p). 



It can thus be seen that these solutions include a step (or shock) gravitational wave if the 

Generally, the Weyl tensor is C n across 



2' 



lowest term in the expansion for V has k 
the wavefront if k m i n = \n + l. 

It can be shown that V and M and their derivatives with respect to rj and p are all 
bounded as p — > for all rj > 0. It thus follows from Q8.5| ) that the component ^2 of the 
Weyl tensor is singular on the axis p = 0. Thus the axis of the cylindrical waves is a scalar 
curvature singularity even though the fluid density on it given by (5.11) is bounded. This 
axis may be considered both as the source of the cylindrical waves, and also as a "sink" 
of the backscattered radiation. 

The possibility of obtaining exact solutions of this type was mentioned in ||] but the 
explicit solutions are given here for the first time. By superposing two sets of terms 
similar to those included in (3J2) but with wavefronts given by 77 = p — a and 77 = b — p 
for positive constants a and b, it is also possible to consider the collision of outgoing and 
incoming cylindrical gravitational waves. This is qualitatively described in figure 3. Exact 
solutions describing this situation are only singular on an initial hypersurface and on the 
axis of symmetry. However, although the qualitative features of such a solution can be 
determined, the complete integral for M cannot be obtained in this case using the methods 
described above. 
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9 Plane waves in an open FRW model 



We now turn to the case in which a plane gravitational wave propagates in an open FRW 
universe with a stiff perfect fluid. In this case, it is convenient to adopt the coordinates 
77, /i, x, y as defined in case 2a of Section 5. In the background region, the metric functions 
are given by and the fluid by Q5.12 ). Since e~ u = smh2r] e~ 2fJ/ in this case, we can 
put 

/ = i(e 2 ^-l), g = \(l-e-^). (9.1) 

We can now consider a gravitational wave with a plane wavefront given by 77 = \i as 
described in Section 3, with the region 77 < /j being part of the exact open FRW model 
into which it propagates. 

Since V = in the background region in this case, the gravitational wave can be 
introduced by putting 

V( V , n) = &( V - M ) Q cf>(k) sinh fc 2r? e~^H k ( ^g 3 " ) dk, (9.2) 

where H k is most conveniently given by ( |6.12| ) and its extension. After applying again the 
method described in ||, the complete integral for M is found to be 

M(t?, fj) = - log sinh 2t? + M G(t? - /j), (9.3) 

where 

rco 1 pn— 1/2 

M = -\\ smh n 2i] e~ 2n ^-dn </>(k)<f>(n - k) 

Ji n J1/2 

~k(n - k)H k H n „ k - ( l - 2c ° s ^f!i + ^ ) fffc-i^*-i] dk, (9.4) 

in which the argument of H k is (e 2/i — cosh 2r?) / sinh 2r/. 

As previously, we can determine the behaviour of the gravitational wave near the 
wavefront where 77 — is small by considering the lowest non-zero term in k in ( |9.2| ). 
Taking this lowest term in the form 

V%/x) = 2* sinh* 2r?e- 2fc ^ fe 6(77 - z), (9.5) 

it can be seen that, near the wavefront, the non-zero Weyl tensor components behave as 

* ! ~-7w('-"» ae ('-"I' ( M > 

4 sinh 27? 
2 sinh'' 2t? 

It can thus be seen that an impulsive gravitational wave is again included here if the 
lowest term in the expansion for V has k = \. It has a step (or shock) gravitational wave 
if the lowest term has k = |. Generally, the Weyl tensor is C n across the wavefront if 

^min — ~\~ 2' 
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From the above results, it can be seen that these waves in open FRW universes are 
qualitatively similar to those in spatially flat FRW universes. They both include im- 
pulsive, shock or smooth fronted wavefronts. As discussed in Section 3, the geometrical 
properties of these plane wavefronts in open FRW universes are described in figure 1 and 
equation (3.7). 

A discussion of these waves and their wavefonts, as well as the possible collision of 
such waves, was given in M. However, the complete solution for M and an analysis of the 
character of the Weyl tensor components near the wavefront has been included here for 
the first time. 



10 Cylindrical waves in an open FRW model 

We now consider a gravitational wave which propagates into a stiff fluid open FRW uni- 
verse along the null cylindrical wavefronts described in figure 2 and Eq. ( 3.1C| ). In this 
case, it is convenient to adopt the coordinates rj, p, (ft, z as defined in case 2b of Section 5. 
In the background region, the metric functions are given by (5.7) and the fluid by ( 5.12| ). 
Since e~ u = sinh2rysinh2p in this case, we can put 



/ = — sinh 2 (ry — p) , g = sinh 2 (77 + p) . 



(10.1) 



We can now consider a gravitational wave with a cylindrical wavefront rj = p as de- 
scribed in Section 3, with the region rj < p being part of the exact open FRW universe 
into which the wave propagates. This can be achieved by putting 



V(r/, p) = logtanh/3 + V(r],p)Q(rj - p), 



(10.2) 



where 



roc / . 

V(V, P) = j i/2 Hk) smh fc 2r, sinh fc 2p H k ( ^X^ 1 ) dk > ^ 

in which H k is now most conveniently given by ( |6~TT|) . In this case the integration of the 
subsidiary equations (^3) is more complicated. However, the complete integral for M can 
still be obtained and expressed in the form 



M(r], p) = — log sinh 2r] + M(ry, p) @(r] — p), 



[10A) 



where 
M 



in which 



1/2 

l 

2 



<t>(k)\2kA k „ l (T,QH k + 



(C 2 -l) 1/2 lr k + (C 2 -l)A k (r,C) 



H, 



k-l 



dk 



00 1 

-r n dn 
1 n 



n-1/2 



1/2 

k(n - k)H k H n _ k 



4>(k)(j)(n — k) 



(C 2 - 



A k (r,C) 



[(2 + a) 2 - * 2 ] 1/2 



n—k—1 



dt. 



dk, 



(10.5) 
(10.6) 
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Hk = -f^fc(C) are as defined above, and r and ( are functions of rj and p given by 



cosh 2rj cosh 2p — 1 

r = sinn 2r/ smh 2p, Q= . (10.7) 

sinh 2rj sinh 2p 

It can be seen that these solutions are singular at the initial big bang at which rj = 0, 
and that there is an additional curvature singularity on the axis of symmetry p = 0. The 
singular axis can be considered to represent both the source of the cylindrical gravitational 
wave and the sink of the backscattered radiation, as in the spatially flat case described in 
Section 8. 

It is also possible to analyse the character of the wave along the wavefront using the 
same method as in previous sections. Accordingly, we include only the lowest non-zero 
term in k in ( 10. 2| ) in the form 



V(V,p) = ^smh k 2i 1 sinh k 2 P H k («£g2g»gi) . (i .8) 

With this it can be shown that, near the wavefront, the non-zero Weyl tensor components 
behave as 

».~- "'^g t " h->)'-^-rt, 

2 smh 2r] 



*2~ T^T (v-p) 2k e(v-P), (io.9) 

2 smh 2rj 



2 smh 2r] 

These expressions can be seen to be very similar to those for cylindrical waves in the 
spatially flat FRW universe as described in Section 8 (see (|8.7| ) in particular). As there, it 
can be seen that these solutions include a step (or shock) gravitational wave if the lowest 
term in the expansion for V has k = n- Generally, the Weyl tensor is C n across the 
wavefront if k m i n = + 1. 

The possibility of obtaining explicit solutions of this type was mentioned in Q but the 
solutions are described in detail here for the first time. It is also possible to consider the 
collision of opposing outgoing and incoming cylindrical gravitational waves just as in the 
spatially flat case. However, although the qualitative features of such a situation can be 
determined, a complete integral for M in the interaction region cannot be obtained in this 
case using the above method. 



11 Toroidal waves in a closed FRW model 

We now consider waves in a closed FRW universe. In this case, no plane- fronted waves can 
exist. However, as described in Section 4, waves with toroidal surfaces can be constructed. 
These are analogous with the cylindrical-fronted waves in an open universe considered 
in the previous section. Mathematically, they can be obtained by replacing the hyper- 
bolic functions by the corresponding trigonometric functions. Physically, the cylindrical 
wavefront becomes closed to form a toroidal wavefront. 
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In this case, it is convenient to adopt the coordinates 77, £, 5, a as defined in case 3 
of Section 5. Starting from the background ( |5.8| ), we may consider the family of null 
hypersurfaces T] — ( = const. At any fixed time r/ = const., the wave surface is a toroid 
C = const. However, as explained in Section 4, two axes now exist — at C = and at 
£ = it/2. We may consider a wave starting at the initial big bang singularity (77 = 0) 
at the circular axis £ = 0, and propagating along the null hypersurface 77 = £ into the 
FRW closed universe. At any fixed time 77, the wavefront will be a toroidal 2-surface as 
illustrated in figure 4. As the universe expands, the front of the wave propagates towards 
the (circular) axis £ = tt/2. 

In the background region, the metric functions are given by (^^) and the fluid by 
fl5.13p . Since e~ u = sin 277 sin 2£ in this case, we can put 

/ = _ sin 2( 7? _ C ) ) 5 = sin 2 (77 + C). (11.1) 

The gravitational wave with the toroidal wavefront 77 = £ is determined by putting 

v( v ,0 = log tan c + v(v,CMv-0, (H-2) 

where 

/OO , v 

J>(k) sin fc 2r, sin fc 2( H k ( l^gggggg ) dfc, (11.3) 

in which iJ^ is given by ( |6.11| ). As above, it can be shown that the complete integral for 
M can be expressed in the form 

M(f],C) = -logsin2r/ + M(?7,C)e(r/-C), (11.4) 

where 

M = -lJ^(k)^2kB k ^ 1 {r,0H k + (f -l) l/2 \r k - (i 2 -l)B k (T,i) H k ^dk 

roo 1 rn— 1/2 

_i / -<>./,-, / 



r n dn / 0(fcWra - k) 
1 n J1/2 



k(n - k)H k H n „ k - {f - l)ff fc _!H n _fc_i dfc, (11.5) 



in which 



and r and £ are functions of 77 and £ given by 

• 1 • 0/ t 1 ~ cos 277 cos 2C , , 

r = sin2r ? sm2C, £ = — — . , (11.7) 

sm 277 sin 2C 

and to avoid confusion with the coordinate £, in this case, we have relabelled the argument 
of H k putting H k = H k (g). 

The character of the wave along the wavefront can also be analysed as in previous 
sections. Accordingly, we consider the lowest non-zero term in k in ( 11.2j ) in the form 

V(V, = ^ sm fc 277 sin* 2 C H k ( ^gffiggt ) . (11.8) 



19 



With this it can be shown that, near the wavefront, the non-zero Weyl tensor components 
behave as 



W4 _ . t (2 fc + 1 )(4 t + l) 
2 sin 2rj 

^ fl fc(2fc + l) {7] _ crQ{7] _ a (1L9) 
2 sin 2r\ 



2 sin 2rj 

These expressions are obvious modifications of those for cylindrical waves in spatially flat 
and open FRW universes as given by ( |8.7| ) and ( 10.9D - They include a step (or shock) 
gravitational wave if the lowest term in the expansion for V has k = i. Generally, the 
Weyl tensor is C n across the wavefront if k m \ n = \n + 1. 

It can be seen that these solutions are singular at the initial big bang at which r/ = 0, 
and at the final big crunch at which r] = tt/2. There is also an additional curvature 
singularity on the axis of symmetry £ = 0, which can be considered to represent both the 
source of the toroidal gravitational wave and the sink of the backscattered radiation. 

As noticed in Section 5, from the big bang until the big crunch, a photon in a closed 
universe with a stiff fluid succeeds in travelling only one quarter of the distance around 
the whole universe. Similarly, an exact shock wave starting at the big bang on the "polar 
axis" C = will reconverge onto the "equatorial" axis £ = tt/2 at the big crunch. The 
space-time diagram of such a situation is given in figure 5. 

It is also possible to consider the collision of two shock waves propagating into the 
closed FRW background with a stiff fluid. At the big bang, one wave starts at the axis 
C = and the other at £ = 7r/2. As the universe expands, the waves approach each other 
and collide just at the time of maximum expansion. One may also consider other initial 
conditions with the waves being some distance from the axes at time 77 = 0. Such waves 
will collide while the universe is still expanding. However, although the qualitative features 
of such a situation can be determined a complete integral for M cannot be obtained in 
this case. Locally, such waves were considered in ||. (The coordinates z,x : y used in |J 
are related to those above by ( = z, 5 = 2V 2 x, a = 2VV) The toroidal character of the 
waves was, however, not noticed and the discussion of the singularities was incorrect. 

Solutions of this type are closed universes with two global, commuting and hypersurface 
orthogonal Killing vectors. They may therefore be considered as Gowdy-type space-times 
|2(|| , pip. It can readily be verified that here the regularity conditions across the null 
hypersurface rj = £ are satisfied. Such conditions were studied carefully by Gowdy himself 
in the vacuum case, and by Kitchingham | |22| ] (see also references therein) in the case of 
stiff fluid cosmologies. 



12 Concluding remarks 

Exact gravitational waves with two isometries have been constructed which can be con- 
sidered as exact perturbations of FRW cosmological models with a stiff fluid. The waves 
have distinct wavefronts and propagate into any of the three possible FRW backgrounds. 
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Although much work has been done on considering various types of waves in the context 
of cosmological models (c.f. reviews [|]]-[||]), waves with distinct wavefronts propagating 
into the standard FRW backgrounds have only recently been considered (||-§§)- 

In the spatially flat and open universes, the wavefronts can be either plane or cylindri- 
cal. In the closed models, they are toroidal. One may consider toroidal waves propagating 
in other models as well, provided that the spatial sections are compactified by making 
appropriate identifications. However, such identifications result in anisotropic background 
spaces. Here we have assumed the standard backgrounds with natural topologies admit- 
ting globally the six-parameter group of isometries. 

We have obtained exact solutions for plane- fronted gravitational waves propagating in a 
spatially flat or open FRW universe with a stiff fluid. These include impulsive waves, shock 
waves, or waves with Weyl tensor components of any arbitrary finite degree of smoothness. 
We have also presented exact solutions for cylindrical- or toroidal-fronted gravitational 
waves propagating in a spatially flat, open or closed FRW stiff fluid universe. These do 
not include impulsive waves, but may include shock waves or waves of any arbitrary finite 
degree of smoothness. 

Solitonic perturbations of FRW models with cylindrical symmetry, without wavefronts, 
have commonly been treated in the literature (see e.g. [14] and 22], and references therein). 
The possibility of constructing plane-fronted waves in the open FRW universes appears 
to have been recognised only recently On the other hand, it is well known that the 
G2 cosmological models, that are often interpreted as inhomogeneous generalisations of 
the FRW models, can sometimes also be used to describe cylindrical perturbations simply 
by reinterpreting the relevant coordinates. It can be seen from the above examples that, 
not only is the topological character of the wavefronts changed by this procedure, but the 
character of the wave may be altered as well. 

No singularities, except those representing the big bang or big crunch, arise in the 
closed vacuum models considered by Gowdy |2(J , |^lj or in the soliton-type waves studied 
by Belinskii [14], Kitchingham [22] and others (see JjJ], 0], 0] for references) in cases of 
both open and closed FRW stiff fluid models. Their explicit solutions, although without 
wavefronts, clearly correspond to our cylindrical waves, or to toroidal waves in the closed 
case. However, in the spatially flat and open FRW models, they could also be adapted to 
infinite plane perturbations. 

With a clear geometrical picture of the wavefronts available, it may be of interest to 
consider also other types of shock waves with two isometries propagating into these FRW 
backgrounds. These may be just test radiative fields (not necessarily gravitational), or 
exact (large) perturbations. For example, one could discuss fluid shocks in all types of 
FRW models, as Tabensky and Taub |l9| considered shock acoustic waves in the spatially 
flat universe in their pioneering work. 

The advantage of studying shock gravitational waves with distinct wavefronts is not 
only in having a clear picture of the wave surfaces. One can also discuss the head-on 
collision of such waves. In the case of the spatially flat and open universes, the collision 
of waves with plane surfaces was analysed in detail in |j] and |7|. The presence of the 
gravitational waves slows down the rate of expansion, but future spacelike singularities do 
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not occur for interacting waves in these expanding universes as they do in the vacuum case 
with a Minkowski background. This result is qualitatively similar to those of Centrella 
|23|1 and Centrella and Matzner |24j which study the collision of waves in an expanding 
vacuum Kasner background. 

Similarly, one can consider the collision of "outgoing" and "incoming" waves with 
cylindrical wavefronts in the spatially flat and open universes, and toroidal wavefronts 
in the closed models. No future spacelike singularities seem to arise again. However, it 
appears that making the wave surfaces finite along one direction in the cylindrical case and 
fully finite in the toroidal case gives rise to timelike singularities representing the histories 
of the corresponding axes. The "source" at the axis appears to be needed to support the 
propagation of such waves. 



Appendix 

We here consider the geometry of possible wave surfaces in a closed FRW universe in 
greater detail. We will show that a typical wave surface is a 2-torus in the 3-sphere 



described by (4.3) and (4.4). First consider a circle S, as illustrated in figure 6, in the 



W, Z-plane given by 

S : W 2 + Z 2 = K 2 Rl, X = Y = 0, (A.l) 

where k € [0, 1] is a constant. The circle S is generally inside the 3-sphere. Next, consider 
a 2-plane given by 

{X, Y)eIR 2 , W,Z fixed, W 2 + Z 2 = n 2 R 2 Q . (A.2) 

Such a plane for any fixed W and Z touches the circle S at just one point, and it intersects 
the 3-sphere in the circle 

X 2 + Y 2 = (I- n 2 )Rl, (A.3) 

with W, Z fixed and satisfying W 2 + Z 2 = k 2 R^. When W = kR and Z = 0, this is the 
circle C illustrated in figure 7. In figure 6 only two points a and b of this circle can be 
seen because one dimension has been suppressed. Now by moving the point (W, X) along 
the circle S and considering the corresponding family of 2-planes (|A.2|) , we find that the 
planes intersect the 3-sphere in the 2-torus T 2 = S 1 <8> S 1 given by 

W 2 + Z 2 = n 2 Rl, X 2 + Y 2 = (1 - k 2 )R,1 (A.4) 

In figure 6, the slice Y = through the 2-torus is seen as two circles, A and B, 



given by (A.4) with Y = 0. These circles represent the intersection of the 2-sphere W 2 + 
Z 2 + X 2 = Rl through the 2-torus flO|). In a universe with large Rq, the circles will 
locally approximate to straight lines and, if k is close to 1, the surfaces will appear to be 
cylindrical. 

Suppressing the coordinate Z, the 2-torus is represented by the circles X 2 + Y 2 = 
(1 — k?)Rq, W = ±kRq, denoted by C and D in figure 7. These two circles represent the 



intersection of the 2-sphere X 2 + Y 2 + W 2 = Rq through the 2-torus flA.4 ). 
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It is seen from (A.4) that the 2-torus is specified by the choice of the parameter k 



which may conveniently be written as n = cos£, with < ( < tt/2. With the form of 



(A.4), the natural coordinates on the 2-torus are the angles, say a and 5, parametrizing 
the two circles of radii Rq cos C and Rq sin£. This leads to the alternative parametrization 
of the whole 3-sphere given by 

W = Ro cos C cos a, 
Z = Rq cos C sin a, 
X = Rq sin C cos 5, 

Y = i? sinCsin5, (A.5) 

where £ € [0, vr/2], a G [0,27r], 5 G [0, 2tt]. In figure 6, the two families of circles A and 

B, each of radius -RocosC have 5 = and 5 = it respectively, and are parametrized by 
< o < 2ir. Similarly, the families of circles C and D, illustrated in figure 7, have radius 
Rq sin£, with a = and a = it respectively, and are parametrized by < 5 < 2ir. Thus, 
this is a natural parametrization for the description of toroidal 2-surfaces inside the 3- 
sphere. Notice that degenerate points occur in both figures either when £ = 0, or when 
£ = tt/2. These are coordinate singularities in which the norms of the Killing vectors d$ or 
d a vanish. It is easily verified that the 3-sphere is covered once with the angular variables 

C, <7, 5 taking values in the given intervals. On the 3-sphere, these angular variables are 
related to the familiar Euler angular coordinates 9, (ft, (which should not be confused 
with the angles used in ((OJ) etc.) by = 2£, <f> = a — 5, ip = a + 5. The Euler angles 
are also often used since the 3-sphere has the same topology as the group SU(2), and the 
generators of the symmetries can be suitably expressed in terms of these angles (see e.g. 
@ or H). 



It is possible to consider the 2-torus ( A.4 ) as dividing the 3-sphere into two solid 



tori. Conversely, S 3 can be constructed by identifying the boundaries of two solid tori. 
This appears to have been first explained in relativity literature by Misner |2(| in his 
analysis of Taub-NUT space, using an inversion. We will now demonstrate this, using the 



parametrization (A.l), thus making apparent the character of the coordinates £, a and 5. 
(See also Appendix B of pQ[ .) 

We first cut the 3-sphere fl4.3| ) into two along the hypersurface Z = 0, or a = and 
a = it, to form two solid 2-spheres. (This is similar to cutting an ordinary 2-sphere along 
an equatorial plane to form two hemispheres which can each be mapped, by stereographic 
projection, onto a disc.) This is illustrated in figure 8, in which the solid 2-sphere on the 
left corresponds to the part of the 3-sphere in which Z > 0, that on the right to Z < 0. 
The two spheres are reconnected by identifying corresponding points on their surfaces. 

The full curves inside the two spheres indicated in figure 8 represent the section of the 
sphere in which 5 = and 5 = it for some fixed value of They can be seen to arise from 
the stereographic projection of the hemispheres with Z > and Z < in figure 6 onto 
the (W, X) plane, the centre of projection being the points W = X = 0, Z = —Rq and 
W = X = 0, Z = +Rq. Introducing coordinates W, X for the projected points of the 
hemisphere with, say Z > 0, W = WRq/(Rq + Z), X = XRq/(Rq + Z), and substituting 
for W, X and Z from (|A75|) with 5 = 0, we obtain W 2 + (X - Ro/sm() 2 = iigcotC- 



23 



For fixed values of (, these are segments of circles, such as A in figure 8, with centres on 
X > 0, bounded by Z = 0. Along these segments a is changing and 5 = 0. As £ — > 7r/2, 
the segments degenerate to the "East pole" on the equator. It can thus be seen that the 
circles A and B in figure 6 are projected into the corresponding circles in figure 8. 

The surfaces of both solid 2-spheres, defined by Z = 0, must be identified in order 
to reconstruct the original 3-sphere. In figure 8, such an identification is illustrated in 
terms of the coordinates C, a and 5. Surfaces inside each sphere on which ( is constant, as 
indicated in figure 8, are joined on the surface of each sphere, dividing the 3-sphere into 
two solid tori. One torus has the closed axis given by ( = and is formed by the axes of 
the two spheres illustrated in figure 8. On this axis a changes from to it in the left sphere 
and continues from tt up to 2ir in the right sphere. Along this axis, the coordinate 8 is 
degenerate and the norm of the Killing vector d$ vanishes. The second torus is evident on 
identifying points near the equators of both spheres. The equators form the second closed 
axis. It is given by £ = ir/2 and is parametrized by < 5 < 2n, while the coordinate a 
is degenerate and the norm of the Killing vector d a vanishes. These two axes, ( = and 
( = 7r/2, are crucial in dealing with toroidal waves in closed FRW universes. 
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Figure 1: The embedding diagram of the section Y = (or = and it) through 
the 3-hyperboloid (3.3) in a 3-dimensional Minkowski space. The 2-planes given by z = 
zo = const., (x,y) G JR 2 (or (X, Y) € IR 2 ) are illustrated as parabolae along which x,y 
vary. 

Figure 2: The embedding diagram of the section Y = (or <p = and it) through 
the 3-hyperboloid (3.3) in a 3-dimensional Minkowski space. The cylindrical 2-surfaces 
given by p = po = const., y £ (— oo,+oo), <\> G [0, 2ir) are illustrated in two hyperbolae 
representing two generators of the cylinder along which the coordinate y varies. The axis 
of the cylinder given by p = is also indicated. 

Figure 3: A space-time diagram illustrating the collision of cylindrical gravitational waves 
in a spatially flat or open FRW universe. Region I is the FRW background in cylindrical 
coordinates (5.5) or (5.7). Region II contains an outgoing cylindrical wave, region III 
contains an incoming cylindrical wave, and region IV is the interaction region following 
their collision. The space-time is only singular on the initial hypersurface i] = 0, and on 
the axis p = 0. 

Figure 4: Toroidal waves propagating into a closed FRW background. The figure shows 
two sections through the toroid as in figures 6 and 7 in the Appendix. The unshaded 
regions represent the exact FRW model, with the wave region illustrated schematically by 
the shaded areas. The singular circular axis C = within the wave region is also indicated. 

Figure 5: This space-time diagram illustrates the propagation of a gravitational wave 
into a closed FRW universe. The coordinates 5 and a are suppressed. The shaded region 
contains the gravitational wave. The null wavefront is toroidal as illustrated in figure 4. 
It starts at a line at the initial singularity r] = and converges to another line at the final 
singularity r\ = it/ 2. There is also a curvature singularity on the initial axis ( = 0. 

Figure 6: A section Y = (<^> = and it) through the 3-sphere of radius Rq showing 
circular cuts A and B through the toroidal surface. The circles S, A and B have radius 
kRcj. 

Figure 7: A section Z = (9 = it/ 2) through the 3-sphere of radius i?o showing circular 
cuts C and D through the toroidal surface. The points a, a' and b, b' are on the circles A 
and B respectively, indicated in figure 6. 

Figure 8: The entire spatial geometry of the closed FRW model at a given time is 
represented by two solid spheres. Points on the surface of upper hemisphere on the left on 
which (7 = must be identified with points on the surface of the upper hemisphere on the 
right on which a = 2tt. Points on the surfaces of both lower hemispheres on which a = it 
must also be identified. Within each sphere, the segments are indicated which arise from 
the section 5 = and it through the spheres. They represent the stereographic projections 
of the circles A and B in figure 6 onto the (W, X) plane. 
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